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The overlap operator clS cl continued fraction* 
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We use a continued fraction expansion of the sign-function in order to obtain a five dimensional formulation of 
the overlap lattice Dirac operator. Within this formulation the inverse of the overlap operator can be calculated 
by a single Krylov space method where nested conjugate gradient procedures are avoided. We show that the 
five dimensional linear system can be made well conditioned using equivalence transformations on the continued 
fractions. This is of significant importance when dynamical overlap fermions are simulated. 



1. Introduction 

Let us start with the explicit form of the mas- 
sive overlap Dirac operator with mass ^, 



^(m) = ^(1 + m) + ^(1 - M)75Sgn(i?(-m)), 



(1) 



where H{—m) denotes some hermitian Dirac op- 
erator with mass — m < of the order of the 
inverse lattice spacing. One method to compute 
the action of the sign-function on a vector is by 
approximating the function by a rational polyno- 
mial 11, 



sgn(a;) « sgn{„^„}(x) 



Pn{x) 



(2) 



where P„ and Qm are irreducible polynomials of 
degree n and m, respectively. The reason for us- 
ing rational approximations is that for a given 
approximation error in the infinity-norm the re- 
quired degree for rational approximations is usu- 
ally much smaller than the one for polynomial 
ones, although both converge exponentially in 
their degree. 

There are various methods to perform such an 
approximation. The methods differ in how they 
interpret the notion of the goodness of the ap- 
proximation. One can use any definition of a 



norm characterising the deviation of two func- 
tions from each other. 

The infinity-norm, for example, is defined as 
the maximal error over a given interval [a,b]. 
There exists an algorithm (Remez algorithm [^) 
to find the best (rational) approximation r(x) to 
a function f{x) with respect to the infinity-norm 
(see also the work of Zolotarjov That is, 

the approximation r{x) to a given function f{x) 
is the rational function of a given degree that 
minimises the maximum value of the quantity 
\f{x) — r{x)\ over the interval [a,b]. 

In order to approximate the sign-function by 
a rational polynomial we choose to find an ap- 
proximation for a constant odd function over the 
interval [e,l]. To do so we approximate on 
the interval [e^, 1] to obtain the rational approx- 
imation r{x) and consider the function xr{x'^). 

The resulting approximation can then be writ- 
ten as 



sgn{2„+i,2m}(a;) = xcq 



(3) 



with real and positive pk and qk ■ This factorised 
form is most stable for numerical evaluation. 

One advantage of using rational polynomials is 
that one can use the partial fraction decomposi- 
tion, e.g. for n = m. 
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Sgn{2n+l,2n}(2;) ^ X I Cq + 



fc=i 



Cfc 



qk 



(4) 
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which allows us to compute the action of the sign- 
function at roughly the cost of a single conjugate 
gradient (CG) inversion of 

2. Continued fractions 

It is well known that rational functions are 
closely related to finite (truncated) continued 
fractions which in turn play an important role 
in the analytic theory of rational functions [^j6|. 
Indeed, we can write any rational approximation 
to the sign-function as a continued fraction |Q 



sgn{2„+i,2n}(a;) = Pox- 



(5) 



(3ix - 



a.2n 



Then we observe that the inverse of eq. (||) is just 
the (1, l)-component of the Schur complement of 
a larger system 



V 



(6) 



Instead of inverting (|4|) using a two-level nested 
CG(-like) procedure we can now find the inverse 
of (|g), i.e. its (l,l)-component, in one go, using 
just one CG(-like) search in a single higher di- 
mensional Krylov space 0]. 

The operation naturally associated with ratio- 
nal functions is the multiplication and inversion 
which, of course, is most trivial. For the overlap 
operator, eq. (^, matters are not so simple since 
the additional term proportional to 75 does not 
commute with sgn(_ff) and its approximations. 
We can, however, still make use of some prop- 
erties of (finite) continued fractions. The most 
important one is contained in the following the- 
orem on equivalence transformations: Two con- 
tinued fractions {ao; {ai, . . . , {a„, Pn}} and 
{ttp; {a'l, I3[}, . . . , {a'„, P'^}} are equivalent if and 
only if there exists a sequence of non-zero con- 
stants c„ with Co = 1 such that 

a'l = c„c„_ia„, n= 1,2,3,..., 

f3'„ = CnPn, n = 0, 1,2, . . . . 



This can most easily be seen by writing down the 
continued fraction after a transformation. 



ClC2Q!2 



(7) 



C2P2 



In fact the coefficients Ci parametrise the equiva- 
lence class of a continued fraction corresponding 
to a given rational polynomial. It is important 
to note that although the equivalence transfor- 
mation leaves the value of the continued fraction 
invariant it affects the spectrum and correspond- 
ingly the condition number of the larger system 
(^). One might therefore hope to find within an 
equivalence class of continued fractions the one 
with the smallest condition number. 

3. Application to the overlap operator 

It is easy to apply the equivalence transfor- 
mation to an overlap operator in which the 
sign-function is approximated by a rational func- 
tion and expanded as a continued fraction. In 
complete analogy to the previous section we 
express the linear system of the hermitian over- 
lap Dirac operator 'y^D{ii)iIj = 6 as a simple 
five dimensional system H^^"' {fij'i' — x where 
* = (i/;,<?!)i,... ,(/.2«+i)^, x= (^0,...,0)'^ are 
five dimensional fermion fields and 77 (/i) is 
a block tridiagonal matrix similar to (^) with 

A+j5,-ClklH, +C2k2 H, . . . , -C2n+lk2n+lH OU 

the diagonal, y/ciAl, ^ciC2, . . . , ^C2„C2„+i with 
trivial Dirac and colour structure on the off di- 
agonals and A± = 1/2(1 ± /i). The ki's are 
uniquely determined by the rational approxima- 
tion of the sign-function while the Ci's are free 
parameters describing the equivalence class of 
continued fractions associated with the given 
rational approximation. 

We observe that an equivalence transformation 
on the continued fraction acts in a similar way 
to block Jacobi preconditioning. Indeed we can 
write c ■ H^^\n;ci = 1) • C where 

C is a block diagonal matrix having the four di- 
mensional unit matrix 1 in the first block and 
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i = I, . . . , 2n+l in the remaining ones. For 
the hnear system under consideration this precon- 
ditioning comes for free since the first component 
of and X are not affected at ah. 

It turns out that within an equivalence class 
the condition number of k{H'-^^), can vary 

by orders of magnitude. An obvious choice for 
the equivalence transformation parameters is Ci = 
1/ki, and indeed, with this choice the matrix is 
already well-conditioned. In order to find the 
optimal values for the Ci in the sense of giving 
the lowest k{H''^^) one can study the matrix in 
the free case where the spectrum can be calcu- 
lated analytically. A non-linear minimisation of 
k{H^^^) then yields a set of optimal values for the 
coefficients q. 

How much is the effective gain in a practical 
application? In order to address this question we 
generated an ensemble of pure gauge field config- 
urations with the Wilson gauge action at /3 = 5.9 
on 6"* lattices. We used the standard Wilson 
Dirac operator Hw{—m) with m = 1.5 as the 
input operator for the overlap and a n = 8 ratio- 
nal approximation for the inverse square root on 
the interval [e^ 1] = [0.00724,1]. This generates 
a five dimensional operator with an extent of 18 
in the fifth dimension. Then we calculated the 
inverse of the hermitian overlap operator ^^D^fi) 
using the five and four dimensional formulation 
with different quark masses /i on a random vec- 
tor. 

Figure |l| shows the average number of Hw- 
matrix multiplications for the inversion in three 
cases: the pole approximation, eq. (^), (4(i over- 
lap), the five dimensional operators with 
Ci = 1/ki (5d matrix) and with the optimal 
Ci's as obtained from the free field approximation 
(5(i matrix optimised). The results indicate that 
it is indeed possible to generate five dimensional 
systems which are well conditioned and that the 
obstacle of having large condition numbers for the 
five dimensional system can be avoided. For small 
quark masses the inversion seems to behave bet- 
ter than the one of the corresponding four dimen- 
sional system. This is of particular importance in 
the context of dynamical overlap fermion simula- 
tions where the inversion of the fermion matrix 
plays a crucial role Another advantage of 
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Figure 1. Number of Hw applications for invert- 
ing the overlap operator in different four and five 
dimensional representations as a function of the 
quark mass /x. 

the five dimensional formulation using continued 
fractions is its obvious suitability for parallel com- 
putations. 
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